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In this article we show that the phase-ordering scaling state for binary fluids is not 
necessarily unique and that local correlations in the initial conditions can be responsible 
for selecting the scaling state. We describe a new scaling state for symmetric volume 
fractions that consists of drops of the one component suspended in a matrix of the 
other. The underlying reason for the existence of the newly observed scaling state 
is that the main coarsening mechanism of binary fluids which is the deformation of 
interfaces by flow is not acting, and this leads to a new scaling law. An initial droplet 
state can be formed by a number of physical phenomena. In a unified description this 
can be undestood as local correlations in the initial conditions. Local correlations with 
length £ are believed to be irrelevant when the typical length scale L of the system 
is large (L > Q. Our result shows that these initial correlations, contrary to current 
thinking, can be important even at late times. 



1.1 Introduction 



Phase-ordering is observed in many spatial systems, including ones that under- 
went spinodal decomposition or mechanical mixing. Examples include a wide 
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variety of systems including magnetic systems, binary alloys, binary fluids, but 
also some reaction diffusion systems]]] and some models of gene evolution^). 

In a seminal paper Hohenberg and Halperin showed that most of these sys- 
tems can be categorized to belong to a small number of categories, and each 
member of a category shows the same universal phase-ordering behavior. In 
particular they distinguished between systems in which the order parameter is 
not conserved, e. g. magnetic systems or genetic systems, and systems in which 
the order parameter is conserved, e.g. binary alloys or binary fluids. Liquid 
systems, such as binary fluids and magnetic fluids, can phase-order not only by 
diffusion but also via flow and therefore they are in a universality class of their 
own. 

The phase-ordering process has been conjectured to observe scaling, i.e. the 
systems are characterized by a single length-scale L{t) such that the morpholo- 
gies of the system scaled but that length-scale are statistically self-similar. The 
time dependence of this length-scale is then known as the scaling law. The 
scaling laws can be deduced by dimensional analysis from the equations of mo- 
tion. For systems without hydrodynamics L(t) scales as t 1 / 2 for systems with 
non-conserved order parameter and as i 1 / 3 for systems with a conserved order 
parameter. 

Because binary fluids have more that one mechanism for domain coarsening 
there is more that one growth law. Off-critical quenches that contain a much 
smaller volume fraction of one component form droplets and this morphology 
grows mostly via diffusion leading to a L(t) ~ i 1 / 3 growth law. Most stud- 
ies, however, deal with critical quenches with equal volume fractions for both 
components. Because of this symmetry bi-continuous morphologies are typically 
formed and both diffusion and hydrodynamic growth mechanisms operate. 

There are, however, exceptions to this general rule. Sometimes dynamical 
asymmetries in the diffusion constant of the two phases or their viscoelastic 
properties can lead to droplet morphologies. A stationary morphology that 
consists of circular droplets has no hydrodynamic pathway of coarsening. Once 
two droplets coalesce they can induce a flow that can induce further coalescence. 
The central idea of this paper is that one of the main hydrodynamics pathways of 
coarsening can be suppressed if the morphology consists of droplets leading to a 
new scaling state. The possibility of the non-uniqueness of the scaling state was 
first suggested by A. Rutenberg, although at the time no persuasive numerical 
evidence could be found |]l2|. 

The droplet morphologies that we require are not unusual and can often seen 
in the early stages of viscoelastic spinodal decomposition [|J |], [Io| as well as in 
systems with an order-parameter dependent mobility. We will use a viscoelastic 
lattice Boltzmann method to create the initial droplet morphology and then use 
a symmetric Newtonian lattice Boltzmann method to evolve the initial droplet 
morphology to verify that this new scaling state exists. 

In viscoelastic phase-separation the droplet morphology consists of low- 
viscosity droplets suspended in the viscoelastic matrix. The opposite morphol- 
ogy is observed in mechanical mixing where the low viscosity material will form 
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the matrix phase. We show that both droplet morphologies are stable under 
phase-ordering, in agreement with experimental results. 



1.2 Numerical method 

For the simulations we use the viscoelastic two-component lattice Boltzmann 
simulation introduced in an earlier paper Briefly, in lattice Boltzmann sim- 
ulations densities that are associated with velocities Vi are streamed on a 
lattice according to the lattice Boltzmann equation 

/i(x + Vi At,i + At) - 

/(x,t) + A^A y - [f°(x, (1-1) 

3 

where ff is an equilibrium distribution, Ay is a collision matrix, and v;Ai is a 
lattice vector. The velocity set for our simulation consists of 17 velocities given 
by {(0,0), (1,0), (0,1), (-1,0), (0,-1), (1,1), (-1,1), (-1,-1), (1,-1), (1,0), 
(0,1), (-1,0), (0,-1), (1,1), (-1,1), (-1,-1), (1,-1)}- Note that the last 8 
velocities are the same as the previous eight velocities. This duplicity allows the 
simulation to have two independent stresses which represent a viscoelastic and a 
purely viscous contribution to the total stress tensor. The two contributions are 
used to produce a Jeffrey's model for the stress (see eqn. fll.SP). The algorithm is 
required to conserve mass and momentum, but not energy. Energy conservation 
is replaced by a condition of constant temperature. The macroscopic density p 
and velocity u are defined as 

i i 

To simulate a two-component mixture we now have to consider the densities of 



the two fluids p± and p2- The first lattice Boltzmann equation (1.1) is now an 
equation for the total density p = p\ + p2 and we introduce a second lattice 
Boltzmann equation for gi to describe the evolution of the density difference 

4> = Pi - P2- 

gi (x + ViAt,t + At) = 

g( x ,t) + — [ s< °(x,t)- 5 i(x ) *)] (1.3) 

T 

where we choose a diagonal collision matrix with a single relaxation time r. 
These densities are only defined on the first nine velocities v^. The density 
difference <f> is given by 

i 

By choosing appropriate equilibrium distributions and an appropriate collision 
matrix, we ensure that the following partial differential equations are being sim- 
ulated up to second order in the derivatives but assuming that the relaxation of 
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the viscoelastic stress a is slow [9 ~ 1/^/e): 

d t p + d x (pu) = (1.5) 

p<9 t u + pu.Vu = -d x P + d x (a v + a) (1.6) 
<J V = ^oo(V(pu) + (V(pu)) T -V.U(5) 

+e 3C V.u<5 (1.7) 

a + 6{(j>)(T {x) = -(^o(^-^oo)(V(pu) + (V(pu)) T ) (1.8) 

9 t + 9 x (0u) = DV 2 » + V.{{<f>/p)V.(P-a)) (1.9) 

where <5 is the identity matrix, a is the viscoelastic stress tensor, cm = 
<9 t (7 + u.Ver — ct.(Vu) — (Vu) t ct is its upper convected derivative, P = 
0.5p + O.OO7(V0V0 - O.5V0.V(M - ^V 2 ^ is the pressure tensor, and /i = 
-0.55^/p + 0.25 ln((p + <p)/(p - </>)) - 0, 007V 2 </> is the chemical potential. The 
parameters is^, £oo, and 9 are determined by the eigenvalues of the collision 
matrix. The values for the parameters were At = 1, r = 1, D = 0.5, £oo = 0.31, 
and i^oo = 0.01. For the low-viscosity phase we used 9 = 0.055, = 0.013 and 
for the viscoelastic phase 9 — 39.5 and v — 1.97. For the symmetric simulations 



of Figure 1.1 we used 9 — 0.055 and Vq m = 0.075. 



1.3 Simulations 

We performed simulations of critical spinodal decomposition of a viscoelastic bi- 
nary mixture in two dimensions where one component is much more viscoelastic 
than the other. These simulations lead to the usual morphologies in which the 
viscoelastic phase is connected and the less viscoelastic phase is dispersed p^j. 
We performed our simulations on a 256 2 and a 1024 2 lattice and after about 
1000 iterations the less viscoelastic phase is completely dispersed, although the 
domains are still highly deformed. We used this morphology as an initial condi- 
tion for a simulation where we make both components purely viscous to examine 
the effects of initial conditions that are not symmetric on symmetric binary fluid 
mixtures. We choose the viscosity such that a system started with a random 
initial condition will be in the viscous scaling state. This allows us to distinguish 
the effect that the morphology created by viscoelastic phase-separation has from 
the effect of viscoelasticity itself in the late-time phase-ordering process. 

In Figure [0] we see a comparison of the phase-ordering from a droplet mor- 
phology (a) and a symmetric initial condition (b). The morphologies for both 
systems are shown after 1000, 2000, 4000, and 8000 iterations. We see that 
the phase-ordering of the morphology with a dispersed phase leads to an even 
more pronounced dispersed phase where almost all domains are circular at late 
times. We see that droplet coalescence occurs frequently with only few domains 
vanishing due to the evaporation-condensation mechanism underlying Oswald 
ripening. The droplet coalescence, however, is not frequent enough to change 
the connectivity of the domains. Instead, we observe that domains become more 
circular on average, suggesting that even for very long times we do not expect a 
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Figure 1.1: A comparison of the phase-ordering of two identical symmetric binary 
mixtures after different early time spinodal decomposition. In (a) a phase-ordering is 
seen from an initial morphology generated by viscoelastic phase-separation in which 
the light component is dispersed (the originally low- viscous component). In (b) we 
see the usual symmetric phase-separation. Fig.(c) shows the length scales L(t) for 
(a) L ~ i 1 / 2 and (b) L ~ t 2 ^ 3 . The apparent growth law in (a) is misleading. A 
more detailed study that examines the growth of the droplet morphology over several 
decades in length shows that the growth-law is L ~ t (see fljll). 
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transition to a bi-continuous morphology. (This is no longer true when we leave 
viscous regime and enter into the inertial regime. Here a return to the inertial 



scaling state is observed ]18|]). In Figure |l-l|(c) we see that the growth law for the 
droplet phase appears to be L{t) ~ t 1 / 2 and is smaller than the L(t) ~ i 2 / 3 seen 



for the symmetric phase-ordering shown in Figure 1.1(b). More recent simula- 
tions of the droplet morphology, however, establish scaling over several orders of 
magnitude and find that the actual growth-law is L ~ t pl|. This emphasizes 
the fact that dynamic scaling analysis that does not cover several decades can 
be misleading. Unfortunately such studies are very computationally demanding 
and exist only in few studies. 

The existence of this second scaling state, distinct from the bi-continuous 
scaling state, is important to understanding the late-time regime of viscoelastic 
phase-separation because, even in the absence of viscoelastic effects, we observe a 
droplet-morphology evolving from the initial morphology created by viscoelastic 
phase-separation. This result is also important for practical reasons in processes 
where late-time morphologies need to be controlled. It is well known that mixing 
of high- viscosity and low- viscosity components by means of mechanical agitation 
leads to morphologies where the high- viscosity phase is dispersed in droplets |1| 
for volume fractions of the low- viscosity component of much less than 50%. This 
effect is enhanced if the high- viscous phase is viscoelastic [0. This leads us to 
consider what the late-time morphology of a phase-ordering system with an early 
morphology created by mechanical mixing would be. 

In order to answer this question we performed a simulation of viscoelastic 
phase-separation and after a droplet morphology had been formed at 1000 time- 
steps, we inverted the properties of the two components. We used this state 
as a model for the morphologies of dispersed droplets of the viscoelastic phase 
in a matrix of the low viscosity Newtonian phase that is typical of mechani- 
cally mixed morphologies. We then continued the simulations and observed the 
phase-ordering behavior of the new morphology. The results of these simula- 



tions are shown in Figure 1.2. In Figure |l.2|(a) the morphology for a viscoelastic 



phase-separation of a 50%-50% mixture is shown. We should emphasize that the 



phase-ordering (see eqn. (l.£)) does not have a composition-dependent diffusion 
constant and therefore no domain shrinkage is observed in these simulations. 
Domain shrinkage can be a very slow process in viscoelastic phase-separation 
that prolongs the spinodal decomposition process and makes it more difficult 
to differentiate the early-stage decomposition and the late-time phase-ordering 
processes. This does not reduce the validity of our results, however, since we 
are only interested in the late-time behavior when the domain shrinkage is com- 
pleted. 



The morphologies shown in Figure 1.2(b) are of a simulation where at 1000 



iterations, after the dispersion was achieved, the viscoelastic and the low viscos- 
ity components were exchanged. We see that the morphology remains a droplet 



morphology, albeit now of the viscoelastic phase. Comparing Figures 1.2 (a) and 
(b) we see that the most important factor in determining the final morphology is 
the early stage dynamics, and that phase-ordering scaling morphologies of both 




Figure 1.2: The phase-ordering of a binary mixture of a viscoelastic phase (dark 
color) and a low viscosity Newtonian phase (light color) is shown in (a) and (b). (a) 
corresponds to a usual viscoelastic phase-separation of a 50%-50% mixture and (b) 
shows the phase-ordering after the viscoelastic and low-viscosity domains have been 
interchanged. Morphologies are shown after 2000, 4000, 8000, and 16000 iterations. 
The scaling behavior of the two mixtures is shown in (c) where for (a) L ~ t ' 2 and 
for (b) L ~ 0.4. As in Fig. [Tl] these power-laws have not been established for several 
decades and should therefore not be taken too literally. 
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dispersed viscoelastic and dispersed low viscosity domains exist. 

From Figure 1.2(c) we see that each scaling state appears to have a different 
growth law. The morphology of dispersed viscoelastic domains grows as L ~ t 0A 
whereas the morphology of dispersed low viscosity domains grows as ~ t 2 after 
near circular droplets have been formed. An anomalous slow growth in vis- 
coelastic phase-separation has first been observed experimentally by TanakaJlTt 
who found L ~ i 015 for a system of high molecular- weight polystyrene/ diethyl 
malonate (4.0 wt.%). It must, however, be emphasized that it is impossible to 
conclusively determine scaling exponents from such a short run. It is conceivable 
that a more extensive scaling analysis could lead to different exponents. 

These simulations also emphasize the difference between a morphology after 
spinodal decomposition and after mechanical mixing. The morphology after 
spinodal decomposition is a dispersed low viscosity phase, whereas the state 
after mechanical mixing has a dispersed viscoelastic phase. Subsequent phase- 
ordering does not change the connectivity of these states, in agreement with the 
conventional wisdom that there is a profound difference in states produced by 
spinodal decomposition and mechanical mixing. 



1.4 Conclusions 

In this article we have shown that more than one scaling state exists for late-time 
spinodal decomposition of two-dimensional binary fluids and that local correla- 
tions in the initial conditions or the early time behavior of a phase-separating 
binary mixture can be very important in selecting one of these scaling states. 
We have also explained the difference between viscoelastic phase-ordering states 
after spinodal decomposition and mechanical mixing. Our results show that 
the volume fraction and the physical properties of a mixture do not select a 
morphology by themselves, but that the morphology of the initial state is of 
paramount importance. This is why viscoelastic phase-separation can lead to 
unusual late-time scaling states even when viscoelasticity is no longer important 
at large length scales. 
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